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ABSTRACT 

We build on recent new evolutionary models of Jupiter and Saturn and here extend our calculations 
to investigate the evolution of extrasolar giant planets of mass 0.15 to 3.0 Mj. Our inhomogeneous 
thermal history models show that the possible phase separation of helium from liquid metallic hydrogen 
in the deep interiors of these planets can lead to luminosities ~ 2 times greater than have been predicted 
by homogeneous models. For our chosen phase diagram this phase separation will begin to affect the 
planets' evolution at ~ 700 Myr for a 0.15 Mj object and ~ 10 Gyr for a 3.0 Mj object. We show 
how phase separation affects the luminosity, effective temperature, radii, and atmospheric helium mass 
fraction as a function of age for planets of various masses, with and without heavy element cores, and 
with and without the effect of modest stellar irradiation. This phase separation process will likely not 
affect giant planets within a few AU of their parent star, as these planets will cool to their equilibrium 
temperatures, determined by stellar heating, before the onset of phase separation. We discuss the 
detectability of these objects and the likelihood that the energy provided by helium phase separation 
can change the timescales for formation and settling of ammonia clouds by several Gyr. We discuss how 
correctly incorporating stellar irradiation into giant planet atmosphere and albedo modeling may lead 
to a consistent evolutionary history for Jupiter and Saturn. 

Subject headings: Stars: Planetary Systems, Planets and Satellites: General, Planets and Satellites: 
Jupiter, Saturn, Equation of State 



1. INTRODUCTION 

Over the past 8 years, nearly 120 giant planets have 
been found in orbit around other stars. These planets 
have added immensely to the regions of parameter space 
in which we find giant planets and we are just beginning 
to understand how these interesting (and often hot!) envi- 
ronments affect the evolution of these objects. (See Hub- 
bard et al. (2002), for a review.) However, as we strive to 
understand these strange new worlds we must remind our- 
selves that our understanding of our closest two gas giants, 
Jupiter and Saturn, is far from complete. As our examples 
of giant planets that will always be the most amenable for 
detailed study, it is of great importance to refine our un- 
derstanding of these planets, and the physics that governs 
them, so that we will have confidence in our understanding 
of more distant giant planets. As radial velocity studies 
reach longer time baselines, we are assured of finding plan- 
ets similar to Jupiter and Saturn, and at similar orbital 
distances. These wider orbital separation planets are also 
more likely to be directly imaged because they are farther 
from the glare of their parent stars. This paper focuses 
on applying recent advances in our understanding of the 
evolution of Jupiter and Saturn to hypothetical extrasolar 
giant planets of various masses and orbital distances. 

Our understanding of the evolution of Jupiter and Sat- 
urn is currently imperfect. The most striking discrepancy 
between theory and reality is Saturn's luminosity. Saturn's 
current luminosity is over 50% greater than one predicts 
using a homogeneous evolution model, with the internally 
isentropic planet radiating over time both its internal en- 
ergy and thermalized solar radiation. This discrepancy 
has long been noted (Pollack et al. 1977; Grossman et al. 
1980; Guillot et al. 1995; Hubbard et al. 1999). Homo- 



geneous evolutionary models of Saturn tend to reach an 
effective temperature of 95.0 K (Saturn's current known 
T e ff) in only ~ 2.0 to 2.7 Gyr, depending on the hydrogen- 
helium equation of state (EOS) and atmosphere models 
used. However, purely homogeneous models appear to 
work well for Jupiter. Figure 1 shows homogeneous evolu- 
tionary models for both planets from Fortney & Hubbard 
(2003) (hereafter Paper I). It has also long been believed 
that the most promising route to resolving this discrep- 
ancy is the possible phase separation of neutral helium 
from liquid metallic hydrogen in the planet's interior, be- 
ginning when Saturn's effective temperature reached ~100 
- 120 K (Stevenson & Salpeter 1977a,b). 

Immiscibilities in two-component systems are common, 
and are the byproduct of the interaction potentials of the 
types of atoms (or molecules) in the mixture. Once the 
temperature of a system becomes low enough, the energy 
of mixing becomes small enough that the Gibbs free energy 
of the system can be minimized if the system separates into 
two distinct phases. One phase contains slightly less solute 
(here, helium) in the solvent (here, liquid metallic hydro- 
gen) than there was initially, and the other phase nearly 
pure solute. This is often termed immiscibility, insolubil- 
ity, or phase separation, and the mixture is said to have 
a miscibility gap or solubility gap. In general, the smaller 
the percentage of atoms in a mixture that is solute, the 
lower the temperature the mixture must attain for immis- 
cibility to occur. A common approximation (see Stevenson 
1979; Pfaffenzeller et al. 1995) for the saturation value of 
x, the number fraction of the solute, is 

x = exp(B - A/k b T), (1) 

where B is a dimensionless constant, ks is Boltzmann's 
constant, T is temperature, and A is a positive, pressure 
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dependent constant with units of energy. As described in 
Paper I B should be close to zero, and A is the increase in 
free energy upon addition of a helium (or whatever atom 
in general) to pure liquid metallic hydrogen. In Paper 
I we showed that whether A increases or decreases with 
pressure has important effects on giant planet evolution- 
ary models. In a solar composition mixture x for helium 
is about 0.085. Oxygen, the 2nd most abundant element, 
is down by a factor of over 100. 

Since helium is relatively abundant in the hydrogen mix- 
ture, the helium (which is predicted to be neutral) will 
perturb the structure of the proton-electron plasma. This 
A constant has been calculated by various methods in the 
papers we will mention below to be ~ 1-2 eV, which for 
solar composition leads to temperatures on the order of 
5000-10000 K for the onset of helium immiscibility. (This 
is also dependent on pressure.) Figure 2 shows in detail 
our current knowledge of the high pressure phase diagram 
of hydrogen and helium (Hubbard et al. 2002). Labeled 
are the current interior adiabats of Jupiter, Saturn, and a 
hypothetical 0.15 Mj planet. Relevant experimental and 
theoretical boundaries are also labeled, as are regions of 
calculated helium immiscibility. 

Salpeter (1973) was the first to note the effects of the 
immiscibility of helium in liquid metallic hydrogen on the 
structure and evolution of a hot, adiabatic, hydrogen- 
helium planet. Stevenson (1975), using perturbation 
theory, performed the first detailed calculation of the 
hydrogen-helium phase diagram, in an effort to map the re- 
gions of pressure-temperature-composition space in which 
helium was likely to become immiscible. His calculations 
roughly agreed with estimates of the current pressures and 
temperatures of liquid metallic hydrogen in Jupiter and 
Saturn's interiors. These calculations indicated that as 
pressure increased, the saturation concentration of helium 
in liquid metallic hydrogen would increase, leading to con- 
stant composition curves that slant down and to the right 
in Figure 2. 

Soon after, Stevenson & Salpeter (1977a, b) performed 
detailed calculations on the dynamics and distribution of 
helium in giant planets. They found that when helium 
becomes immiscible in liquid metallic hydrogen, the com- 
position that separates out is essentially pure helium, and 
this helium on fairly short timescales (relative to the con- 
vective timescale) will coalesce to form helium droplets. 
These droplets, once they reach a size of ~ 1 cm, will at- 
tain a Stokes velocity greater than the convective velocity 
and will then fall down through the planet's gravitational 
field. If the droplets reach a region where helium is again 
miscible at higher concentration, they will redissolve, en- 
riching the deeper regions of the planet in helium. They 
found that this "helium rain" could be a substantial ad- 
ditional energy source for giant planets. Helium would be 
lost from all regions with pressures lower than the pres- 
sures in the immiscibility region, since the planet is fully 
convective (or nearly so) up to the visible atmosphere. Ex- 
cess helium would be mixed down to the immiscibility re- 
gion and be lost to deeper layers. This would leave all 
molecular regions up to the visible atmosphere depleted in 
helium. 

Relatively few studies have been done since then on 
phase diagrams of hydrogen-helium mixtures. Hubbard 



& Dewitt (1985), using a Monte Carlo technique, but sim- 
ilar assumptions to that of Stevenson (1975), obtained es- 
sentially the same results. The most recent calculations 
(Pfaffcnzcllcr ct al. 1995) utilized molecular dynamics to 
predict a helium-immiscibility region with a shape very dif- 
ferent from that of Stevenson (1975) and Hubbard & De- 
witt (1985). They find that as pressure increases, the sat- 
uration concentration of helium in liquid metallic hydro- 
gen decreases, which leads to constant composition lines 
that slant up and to the right in Figure 2. Interestingly, 
these constant-composition lines run nearly parallel to the 
giant planet adiabats. Detailed inhomogeneous evolution- 
ary models including helium phase separation were not 
performed until Hubbard et al. (1999). These authors in- 
vestigated the cooling of Jupiter and Saturn when the mass 
of helium rained out linearly with time since Jupiter and 
Saturn's formation, or alternatively, rained out just before 
the planets reached their known effective temperatures. 
These are two logical limiting cases. 

There is observational evidence that helium phase sep- 
aration has begun in both Jupiter and Saturn. The pro- 
tosolar mass fraction of helium is calculated to be near 
F=0.27 (Loddcrs (2003) puts the number at 0.2741). The 
atmospheres of both Jupiter and Saturn are depleted rela- 
tive to this value. With the assumption that these planets 
globally contain the protosolar Y, the missing helium must 
be in deeper layers of the planet. The case for Jupiter's de- 
pletion is clear cut, with a value of Y = 0.234 ± .005 from 
the Helium Abundance Detector (HAD) on the Galileo 
Entry Probe (von Zahn et al. 1998). The case for Sat- 
urn is much less clear. Without a past or planned Sat- 
urn entry probe, Saturn's atmospheric helium abundance 
can only be obtained through indirect methods, using in- 
frared spectra with or without radio occultation derived 
temperature-pressure (T-P) profiles. Analysis of Voyager 
measurements indicated Y = 0.06 ± 0.05 in Saturn (Con- 
rath et al. 1984). However, the mismatch between the 
Voyager derived value for Jupiter (Y = 0.18 ±0.04, (Gau- 
ticr ct al. 1981)) and the accurate HAD measurements, 
along with Hubbard et al. (1999) evolutionary and Guil- 
lot (1999) static models, led Conrath & Gautier (2000) to 
perform a reanalysis of the Voyager data. The details of 
their investigation will not be described here, but by dis- 
regarding the occultation derived T-P profile, which may 
be in error, they obtain Y = 0.18 — 0.25 for Saturn's at- 
mosphere. 

Noting the clear need to better understand Jupiter and 
Saturn in light of these atmospheric Y values, in Paper I we 
calculated the first evolutionary models that coupled high- 
pressure phase diagrams of hydrogen-helium mixtures and 
a grid of radiative atmosphere models for giant planets. 
A variety of Saturn evolutionary models were calculated 
that included helium phase separation. The main findings 
of Paper I were as follows. The phase diagram of Hubbard 
& Dewitt (1985), which is essentially the same as that of 
Stevenson (1975), is inapplicable to the interiors of Jupiter 
and Saturn, if helium phase separation is Saturn's only ad- 
ditional energy source. These phase diagrams predict that 
A from equation (1) is a decreasing function of pressure. 
As Figure 3 shows, this phase diagram prolongs Saturn's 
cooling 0.8 Gyr, even in the most favorable circumstance 
that all energy liberated is available to be radiated, and 
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does not instead go into heating the planet's deep interior. 
Paper I found that if one were to match the Y atmos of Con- 
rath & Gautier (2000) and prolong Saturn's evolution to 
a T e ff of 95.0 K at 4.56 Gyr, the helium that becomes 
immiscible and rains down to deeper layers needs to rain 
far down into the planet, likely all the way to the core, in 
order for enough energy to be released and still match the 
relatively high (Y" atmos =0.18 - 0.25) helium abundance. 

In Paper I an ad-hoc phase diagram was created that 
was essentially a modification of the phase diagram cal- 
culated by Pfaffcnzcller et al. (1995). (Pfaffenzeller et al. 
(1995) find a phase diagram in which A from equation 
(1) is an increasing function of pressure.) To simplify the 
evolution of the planets the ad-hoc phase diagram was 
constructed such that helium immiscibility region runs ex- 
actly parallel to the planets' internal adiabats. Therefore, 
there is no region for the helium droplets to redissolve 
in the liquid metallic hydrogen. This causes helium that 
phase separates to rain all the way down to the planet's 
heavy element core. Consequently, all hydrogen/helium 
regions, molecular and metallic, become more helium poor 
as the helium layer on top of the core grows. This ad-hoc 
phase diagram allows Saturn to reach an age of 4.56 Gyr 
and T e ff 95 K while its Y atm os drops to 0.185. Figure 3 
shows the evolution of T e ff vs. time. With this phase di- 
agram, Jupiter evolves homogeneously to the present day 
and reaches ^4.7 Gyr at 124.4 K without helium becoming 
immiscible. Jupiter would then begin to evolve inhomoge- 
neously at T e f / below 123 K. Still in need of explanation is 
Jupiter's depletion of helium (and neon, which may have 
been carried away in the helium (Roulston & Stevenson 
1995)) in its atmosphere. 

What we have from Paper I is a high-pressure hydrogen- 
helium phase diagram that is calibrated to Jupiter and 
Saturn. Specifically, both planets reach their known ef- 
fective temperatures after ~ 4.6 Gyr, with an improved 
but still imperfect understanding of helium phase separa- 
tion. The purpose of the present paper is to investigate 
the effects of helium phase separation on the evolution of 
theoretical giant planets in orbits around other stars using 
the phase diagram derived in Paper I. 

Here we concentrate on planets ranging in mass from 
0.15 Mj (half Saturn's mass) to 3.0 Mj, and derive the lu- 
minosity, T e ff, radius, and Y atmos as a function of time 
during the planets' evolution. Of the many extrasolar 
planets found to date, the planet with the smallest mini- 
mum mass is HD 49674b, at 0.12 Mj (Butler et al. 2003). 
As we will discuss later, a planet with a mass this small 
likely contains no liquid metallic hydrogen, only dense 
molecular hydrogen. Approximately 70% of all known 
planetary candidates have minimum masses of less than 
3.0 Mj, so planets of the masses we explore here are sure 
to be abundant. 

Our standard models (to be discussed in Section 3.1) 
incorporate a primordial 10 Me heavy element core for all 
planets, but later in Section 3.3 we investigate the effects of 
20 Me and coreless models. (1 Jupiter mass, Mj, is 317.7 
Earth masses, Me-) We will later find that for the more 
massive planets, varying the core mass has little effect on 
T e f f but a large effect on planetary radii. In Section 3.2 we 
investigate the effects of modest stellar irradiation. Stellar 
heating retards a planet's cooling, and we find that if a 
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planet is within a few AU of its parent star, the planet 
will reach its equilibrium temperature before its interior 
reaches temperatures cool enough for helium to become 
immiscible. We calculate the cooling of planets in isola- 
tion and at 10 and 5 AU from a constant luminosity 1.0 Lq 
star. In Section 4 we discuss the atmospheric properties 
and detectability of these EGPs likely to be undergoing 
helium phase separation. 

2. DESCRIPTION OF THE MODELS 

Our evolutionary models are the same as those described 
in Paper I and use the method first described by Hubbard 
(1977). A thermal history calculation for an isolated non- 
rotating giant planet of mass M, radius a, and specified 
composition X (where X is a matrix of the various mass 
fractions of elements), yields relations of the form 

L = 47rcra 2 T e 4 // = L(M, t, X), (2) 

a = a{M,t,X), (3) 

where L is the planet's luminosity, a is the Stefan- 
Boltzmann constant, T e // is the planet's effective temper- 
ature, and t is the planet's age (i.e., time since accretion 
of its hydrogen envelope). Under the assumption of ho- 
mogeneous evolution, i.e., that X(r) = constant, (where r 
is the radius of a mass shell inside the planet) and there- 
fore S(r) = constant (where S is the specific entropy of 
the deep interior), expressions (2) and (3) can be derived 
with the help of a grid of model atmospheres. The grid 
is obtained by choosing independent variables T e // and 
g (the atmosphere's surface gravity), integrating the at- 
mospheric structure inward to a depth where it is fully 
convective and essentially isentropic, and then calculating 
S of the atmosphere at depth: 

S = S(T eff ,g,X), (4) 
where the surface gravity is given by 

g = GM/a 2 . (5) 
For the hydrogen-helium envelope of the planet we use 
the "interpolated" equation of state (EOS) of Saumon 
et al. (1995). Our initial helium mass fraction is Y — 0.27. 
For the heavy element core we use the ANEOS olivine 
EOS. For the heavy elements in the envelope (we take 
Z = 0.02) we use the ANEOS water EOS. We use the 
model atmosphere grid of Burrows et al. (1997), as up- 
dated for Hubbard et al. (1999). It is a non-gray grid 
of atmospheres at solar composition, suitable for isolated 
giant planets and brown dwarfs. We couple our phase dia- 
gram from Paper I to our evolution models in the following 
way. Using our phase diagram and a constant reference 
pressure of 2 Mbar, we calculate the specific entropies of 
adiabats for many Y values, from 0.27 to 0.055, in steps of 
0.005. Our helium poor adiabats then have homogeneous 
helium-poor composition from the outer boundary (1 bar) 
until the beginning of the pure helium region. The number 
of mass shells needed for the pure helium region is calcu- 
lated under the constraint that helium must be conserved 
to 0.05% or better for the entire planet. The entropy of 
the pure helium region is set by the constraint that there 
be no discontinuity in temperature across the boundary. 
The mass of the heavy element core remains the same for 
all models. 

We make several approximations, most of which will af- 
fect our results to fairly minor degrees. The first is that the 
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planets are spherically symmetric and nonrotating. This 
is entirely justified as these planets are hypothetical and 
the rotation rates of non-tidally locked EGPs are entirely 
unknown. For Jupiter and Saturn models, neglecting ro- 
tation introduces errors in evolution timescales of only a 
few percent. We also assume that the primordial heavy 
element cores do not take part in the planets' evolution. 
The vast majority of a giant planet's thermal reservoir is 
in the liquid metallic hydrogen, which has a high heat ca- 
pacity. The energy generated by radioactive decay in the 
core contributes negligibly to the heat budget of a giant 
planet (Hubbard 1980). 

Our treatment of irradiation from the parent star is only 
approximate. Recently, a move has been made to consis- 
tently incorporate the effects of stellar irradiation on the 
evolution of giant planets. This has been shown to be of 
critical importance for close-in giant planets such as 51 Peg 
b and HD 209458b (see Guillot & Showman 2002; Baraffe 
et al. 2003; Burrows et al. 2003). For the planets we study 
here, which are at least 5 AU from their parent star, we 
incorporate the effects of irradiation using an approxima- 
tion that has commonly been used for Jupiter and Saturn. 
We assume the incoming stellar radiation that is scattered 
does not change the planetary atmosphere's T-P profile 
from that of an isolated object. Furthermore, the stellar 
radiation that is absorbed (thcrmalized) is absorbed down 
in the convective region of the planet's atmosphere, which 
shares the adiabat of the deep interior. In this way we 
can use a grid of model atmospheres for isolated objects. 
Following Hubbard (1977), equation (2) must be modified 
as below. 

L ee 4iraa 2 (T* ff - T* q ) = L{M, t, X), (6) 

where T eq is the effective temperature that the planet 
would have if it had no intrinsic luminosity (L = 0). We 
derive T eq from the Bond albedo, As, according to 

^aa 2 T^ q = (1 - A B )7ra 2 £*/47rd 2 , (7) 

where is the stellar luminosity and d is the star-planet 
distance. Using equation (6), an equation can be written 
that gives the heat extracted from the planet's interior per 
unit time: 

f 1 f)Q 

L(M,t,X) = -M dmT—, (8) 

JO &t 

where the dimensionless mass shell variable m is defined 

by 

1 f T 2 

to = — / 4irr' 2 dr'p(r'). (9) 

Put simply, a planet's luminosity is derived directly from 
the decrease in entropy of its interior. Equation (8) is 
valid for cither homogeneous or inhomogeneous evolution, 
because S is a function of the mass shell. (For homoge- 
neous evolution, since X does not vary with the mass shell, 
neither does S.) Regions of the interior that are helium 
rich have a lower entropy per unit mass. As the planet 
cools and the helium region on top of the core grows, 
more mess shells are "converted" from hydrogen-helium 
mixtures (high entropy) to pure helium (low entropy). We 
can then calculate the time step between successive models 
in an evolutionary sequence. Equation (8) can be rewrit- 
ten as 

M f 1 

8t=-— / TdSdm, (10) 
L Jo 



where dt is the time step, and the other variables have the 
same meanings as defined earlier. 

In our calculations we assume that there is no disconti- 
nuity in temperature across the boundary between the he- 
lium depleted and pure helium regions. From the Saumon 
et al. (1995) EOS we calculate the adiabatic temperature 
gradient in the pure helium region and ensure that the 
temperatures of the two regions match at the composition 
boundary. In general this prescription may be too simple. 
As discussed in detail in Stevenson & Salpeter (1977b) and 
also in Paper I, a helium composition gradient will lead to 
an increased temperature gradient in order to maintain 
convective instability. Here we assume an infinitely thin 
boundary between our two regions, but that is an artifact 
of the shape of our phase diagram. More generally, for 
arbitrary phase diagrams, helium composition gradients 
would occur with effects on the temperature gradient. 

3. RESULTS OF CALCULATIONS 

3.1. Evolution of Isolated EGPs 

The evolution of the low-mass EGPs we explore here 
can be drastically changed by the onset of phase sepa- 
ration of helium from hydrogen. Figures 4, 5, 6, and 7 
taken together map out our calculations for the luminosity, 
T e ff, Y a tmos, and radii of isolated EGPs. Solid lines show 
inhomogeneous evolutionary models incorporating helium 
phase separation, while dotted lines are for homogeneous 
models. A planet's luminosity falls at a much smaller rate, 
the gradual contraction of the planet is slowed, and the 
Y a tmos falls as helium is lost from molecular and metal- 
lic regions, with the excess helium raining down to the 
planet's core. 

The lowest mass planet we model is 0.15 Mj, which is 
47.7 Me- This means ~75% of the planet's mass is hy- 
drogen/helium envelope. This is approximately the lowest 
mass a planet can be while still having some liquid metal- 
lic hydrogen in its deep interior. For comparison Uranus 
and Neptune are 14.5 and 17.1 Me, respectively. This 
lower mass limit is dependent on the mass of the core and 
the exact pressure(s) at which hydrogen becomes metallic. 
Later on we will show that coreless models for planets of 
this mass do not contain liquid metallic hydrogen. For our 
standard 0.15 Mj model, the onset of helium separation 
is fairly early in the planet's evolution, at an age of 700 
million years. Of course our determination of this age for 
the onset of phase separation is dependent on the temper- 
atures predicted by the Saumon et al. (1995) EOS and the 
planets' assumed initial helium mass fraction. 

The evolution of planets of increasing masses proceed 
similarly. Note that except for Figure 4, all figures show 
time linearly in Gyr. This makes it clear how quickly giant 
planets fall from their initial large, highly luminous state. 
At a given age, the higher the mass of the planet, the 
higher the entropy of the planet's interior. Consequently, 
the greater the planet's mass, the later in the planet's evo- 
lution helium phase separation begins. For our phase dia- 
gram, the liquid metallic region of the planet must reach 
a specific entropy of 6.11 fcs/baryon before helium begins 
to separate. 

Figure 4 shows the full evolution, starting from 1 Myr, 
of our range of planets. Once helium phase separation 
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begins, the luminosity of the lower mass planets (0.15 to 
0.3 Mj) is increased by a maximum factor of 2.25 over 
the prediction of the homogeneous models. For the higher 
mass planets this ratio is near 1.7. The explanation for 
the effect being larger for the smaller mass objects can be 
understood coarsely in the following way. For each planet, 
the same percentage of the planet's mass is falling onto 
the core. A massive planet (say 2.0 Mj) that is 10 times 
more massive than a 0.2 Mj planet is only ~ 30% larger 
in radius. This means that the helium is falling a compa- 
rable distance, even though the masses differ by a factor 
of ten. The planetary radius increases much more slowly 
with planet mass than the actual physical mass of helium 
raining down. Consequently, the lower mass planets are 
affected to a larger degree. The cooling rate in T e ff per 
Gyr decreases by a factor of 4.5 upon the onset of phase 
separation. A 1.0 Mj planet cools at a rate of -12 K Gyr -1 
just before the onset of phase separation and -2.6 K Gyr -1 
after. This new cooling rate is maintained to within ~ 10% 
past an age of 10 Gyr. 

Changes in planetary T e ff can be on the order of 10-15 
K, when compared to homogeneous models, which is quite 
a marked difference, while radius increases are on the or- 
der of 1000-2000 km, which is somewhat small. However, 
for the 1.5 and 2.0 Mj objects the onset of phase sepa- 
ration leads to a near halt of contraction of the planets. 
Figure 6, which shows Y atm os as a function of time gives 
one a feeling for how helium phase separation is proceed- 
ing in the planet as a function of time. Unfortunately, 
Yatmos will never be an observable quantity except for in 
Jupiter and Saturn. Notice that the 1.0 Mj model begins 
to lose helium to deeper layers at an age of 3.25 Gyr, 1.3 
Gyr earlier than the age of the solar system. The onset of 
phase separation will be delayed for all planets once the 
absorption of stellar photons is accounted for. 

3.2. Effects of Stellar Irradiation 

The proceeding discussion is only correct for planets in 
isolation, such that no energy is absorbed from a parent 
star. Very low mass planets in isolation may exist, but 
Gyr-old isolated planets seem unlikely targets for direct 
detection and characterization. In our solar system, both 
Jupiter and Saturn reradiate more absorbed solar energy 
than internal energy. 

In order to investigate how stellar irradiation may alter 
the effects of helium phase separation, we place our same 
model planets at a distance of 10 and 5 AU from a star 
with a constant luminosity of 1.0 L Q . We make use of 
theory of Hubbard (1977) outlined in Section 2. The ac- 
curacy of this method breaks down when a planet's T e // 
reaches its T eq , however. As has been discussed in rela- 
tion to the intensely irradiated EGPs such as HD 209458b 
(Guillot & Showman 2002), the Hubbard (1977) theory, 
when making use of a model atmosphere grid for isolated 
planets, predicts a halting of the contraction of a planet at 
T e f f = T eq (and hence a halting of the cooling of their in- 
teriors) but in reality a planet likely continues to contract 
while maintaining a constant T e // = T eq , as the radiative 
region of the planet's atmosphere expands to encompass 
higher pressures. This should be kept in mind when inter- 
preting Figures 9 and 11, which show a likely unphysical 
halting of phase separation. The calculations of T e ff in 



Figures 8 and 10 are accurate, as T e ff will eventually fall 
to the T eq value. 

Currently, Jupiter and Saturn have essentially the same 
bond albedo of 0.343. To date, evolutionary models of 
these planets have used this value for their entire evolution, 
which is overly simplified. Since wc are not calculating 
consistent atmosphere models for various stellar distances 
and planet gravities, we too assume our model planets 
have a bond albedo of 0.343. Although progress in mod- 
eling giant planet albedos has been made (Sudarsky et al. 
2000), great uncertainties in the modeling of equilibrium 
condensate clouds in EGP atmospheres remain (Ackerman 
& Mar ley 2001; Cooper et al. 2003). Since an understand- 
ing of nonequilibrium condensates is needed to account for 
the albedos of Jupiter and Saturn, much work still needs 
to be done. Prospects for an accurate understanding of 
EGP albedos in general will be helped by detections of 
reflected light from EGPs. 

Since the store of internal energy for the lowest mass 
planets is the smallest, and they consequently have the 
lowest T e ffS at a given age, they will be the most affected 
by stellar irradiation. At a distance of 10 AU, a planet's 
T eq is 79.2 K and at 5 AU it is 112.0 K. For reference, 
if one were instead to consider a M5V star of luminosity 
L=0.22 L Q , these T eq s would correspond to distances of 
2.3 and 4.6 AU. 

Figures 8 and 9 show how the planets' T e ff and Y atmos 
are affected by stellar irradiation at 10 AU. The onset of 
helium phase separation is delayed by ~ 300 Myr, but the 
planetary evolution proceeds in a fashion similar to the 
isolated case for the more massive planets. However, the 
less massive planets are qualitatively affected by the stel- 
lar heating. The 0.15 and 0.2 Mj planets reach their T eq 
in 7 and 8 Gyr, respectively, and at that point their energy 
budget is dominated by the stellar, rather than intrinsic, 
source. The age of the 1.0 Mj planet at the onset of he- 
lium phase separation is 3.5 Gyr, still only ~ 77% of the 
age of the solar system. 

The effects of stellar irradiation are even more clear 
in Figures 10 and 11. These figures show planets with 
T eg =112.0 K. The planets with masses less than 0.7 Mj 
reach their T eq before helium phase separation can even 
begin. With their evolution stalled at this early age, our 
theory predicts that they never cool enough for helium to 
become immiscible in their interiors. However, as these 
planets likely do continue to contract at T e ff = T eqi a 
more correct statement would be that they begin phase 
separation well after their T e f / is dominated be absorbed 
stellar energy rather than internal energy. The time of on- 
set of phase separation is extended about 2 Gyr for all the 
planets. The 0.7 Mj planet has just begun to lose helium 
to deeper layers when its evolution is stalled, as is seen in 
Figure 11. For the 1.0 Mj object, helium phase separation 
begins at an age of ~ 5.4 Gyr. 

For both sets of irradiated models, the planetary radii 
can be calculated by reading off the mass and T e f f from 
Figure 8 or 10, consulting Figure 5 for the age of an iso- 
lated model for that M/T e ff, and then reading the radius 
off Figure 7 at the isolated planet age. This works because 
of our approximation that the stellar flux is absorbed in 
the convective region of the planet. 
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3.3. Effects of Alternate Core Masses 

The current mass of Jupiter's and Saturn's heavy el- 
ement cores can be constrained with static models that 
must match each planet's mass, rotation rate, radius at 
1 bar, temperature at 1 bar, and gravitational moments 
J2, J4, and Jg, using given EOSs for hydrogen, helium, 
and heavier elements. The most recent estimates for these 
core masses, taking into account uncertainties in all pa- 
rameters, give a core mass of 0-10 Me for Jupiter and 
10-20 M E for Saturn (Saumon & Guillot 2003). 

Our understanding of Jupiter and Saturn led us to 
choose 10 Me as a realistic core mass for our standard 
models. Prospects for determining whether or not EGPs 
possess cores are uncertain. For transiting EGPs, if a de- 
tection of a planet's oblateness could be measured from a 
light curve, and its rotation rate derived with some other 
method, some indication of the amount of central concen- 
tration of the mass could be estimated (Barnes & Fort- 
ney 2003; Seager & Hui 2002). A more promising route 
may come from obtaining the radii of low-mass transit- 
ing giant planets, where differences in radii for planets 
with and without cores will be large (Bodenhcimer et al. 
2003). However, a number of unknowns may complicate 
this picture, as uncertainties in stellar age, radius, and 
mass translate directly into uncertainties in planetary age, 
radius, and mass. These effects may be somewhat difficult 
to untangle. 

Therefore, it is worthwhile to investigate what effects 
larger or smaller core masses may have on the evolution 
of EGPs. Our alternate models have core masses of zero 
and 20 Me- As one would presume, and can be seen in 
Figure 12, the most obvious effect is on the radii of the 
planets. A larger mass of heavy elements will significantly 
reduce a planet's radius. This difference can be 10,000 
km for the lower mass planets. (This is true whether the 
heavy elements are in the core or are uniformly mixed.) As 
the model planets increase in mass, the difference between 
zero and 20 Me of heavy elements becomes a decreasing 
smaller percentage of the planets' mass, and therefore has 
a corresponding smaller effect on the planetary radii. 

As was noted in Figure 2, EGPs below a certain mass 
will not have a high enough central pressure to possess any 
liquid metallic hydrogen. Therefore, these planets can- 
not undergo phase separation of helium at these temper- 
atures, because no liquid metallic hydrogen exists for the 
helium to become immiscible in. (Helium likely becomes 
immiscible in molecular hydrogen at some temperature, 
but this is almost assuredly at temperatures much lower 
than those found in giant planets (Stevenson & Salpeter 
1977a), except perhaps near ~ 1.0 Mbar.) At a given 
planet mass, the greater the mass of the core the larger the 
central pressure in the core and the larger the pressure at 
the core/envelope boundary. Since the pressure at which 
hydrogen turns metallic is currently unknown we simply 
choose 2 Mbar as the transition pressure, independent of 
temperature. Our coreless 0.15 and 0.2 Mj planets do not 
reach this central pressure in 10 Gyr, while the coreless 0.3 
Mj planet does in ~ 1 Gyr, so helium immiscibility pro- 
ceeds as described earlier. (Current evidence indicates the 
hydrogen insulator/metallic transition is likely continuous, 
rather than 1st order, and theory predicts this transition 
is also a function of temperature. See Figure 2 and refer- 



ences in the caption.) 

This effect leads to an interesting bifurcation in the evo- 
lution of the lowest mass EGPs. The 0.15 and 0.2 Mj 
coreless planets possess no liquid metallic hydrogen and 
cool homogeneously during their entire evolution, unaf- 
fected by helium phase separation. The models with 10 
and 20 Me cores are affected by helium phase separation, 
and consequently have quite different T e //s after several 
Gyr. This effect can be seen in Figure 13 which shows the 
evolution of coreless EGPs and those with 20 Me cores. 
The models with large cores can have T e ffS ~ 15 K higher 
than coreless models — a difference of ~ 20% due to core 
size alone! The coreless 0.15 and 0.2 Mj planets are shown 
as dashed lines for clarity. It should be noted that if he- 
lium immiscibility docs occur at planetary temperatures 
in very dense conducting molecular hydrogen near ~ 0.7-1 
Mbar, this bifurcation in evolution would occur at slightly 
lower masses. 

4. DISCUSSION 

4.1. Ammonia Cloud Formation Timescales 

From the proceeding discussion and figures it is clear 
that helium phase separation can be a substantial addi- 
tional energy source in EGPs, therefore making these ob- 
jects more luminous than one would predict from homo- 
geneous models. This increased luminosity will also delay 
the time until the formation of condensates in the planets' 
atmospheres. For our phase diagram considered here, since 
helium phase separation does not begin until the EGPs are 
rather cold, most abundant condensates will have formed 
and moved to high pressures (well below the visible at- 
mosphere) before helium becomes immiscible in a planet's 
interior. However, the formation of ammonia clouds is 
during or after the onset of helium phase separation, and 
consequently, timescales associated with these clouds will 
be delayed due to this additional energy source. This will 
change the time these planets are Class I EGPs, as defined 
by Sudarsky et al. (2000) and Sudarsky et al. (2003). 

We calculated the time necessary for our planetary adi- 
abats, at a pressure of 1 bar, to reach the condensation 
curve of ammonia, as taken from Ackerman & Marley 
(2001). This is the time necessary for the ammonia cloud 
base to reach 1 bar. We were interested in whether the 
additional energy source due to helium phase separation 
could delay or stall the formation of clouds in a planet's 
atmosphere. This calculation is only correct if the planets 
are adiabatic by a pressure of 1 bar in the atmosphere. 
This seems to hold true for Jupiter and most EGP atmo- 
sphere models for isolated or moderately irradiated plan- 
ets (Sudarsky et al. 2003). Since we do not compute at- 
mosphere models, we cannot describe the effects on cloud 
formation at higher altitude, lower pressure regions of the 
atmosphere that do not lie on the interior adiabat. As the 
pressure decreases, atmospheric T-P profiles become more 
isothermal. But here we outline the general effect. 

Figure 14 shows, for our isolated and irradiated models, 
how phase separation and stellar irradiation change the 
time for the ammonia clouds to reach 1 bar. The delay 
due to phase separation is always greater than ~ 1 Gyr, 
and can be 6 Gyr or longer for more massive planets. Plan- 
ets within 5 AU of the parent star only marginally reach 
this point with no helium phase separation. The lowest 



Phase Separation in Giant Planets 



7 



mass planets reach their T eq before the ammonia clouds 
can reach 1 bar. When phase separation is included this 
is true for all planets, independent of mass. What re- 
mains unclear without atmosphere models is whether he- 
lium phase separation delays the start of the formation of 
the ammonia clouds, or alternatively, if the clouds have al- 
ready begun forming before phase separation begins, and 
the helium phase separation merely greatly extends the 
timescale the clouds are in the visible atmosphere. If we 
use Saturn as a guide, which has visible clouds that are 
thought to be composed of condensed ammonia, it is likely 
the latter, although the position of the condensation curve 
of ammonia will depend on the planet's nitrogen abun- 
dance. It seems that due to the predicted delay in cooling, 
ammonia clouds will reside in a planet's visible atmosphere 
for several additional Gyr. 

4.2. Applications to the Cooling of Jupiter and Saturn 

As can be seen from the models in Section 3, the amount 
of energy absorbed by a giant planet due to stellar irradia- 
tion has a critical effect on the time for the onset of helium 
phase separation. The model planets at 5 AU reach helium 
phase separation 2 Gyr later than the isolated planets. 
The fact that we highlighted the evolution of 1.0 Mj plan- 
ets especially is no accident. Jupiter's Y atm os is only 86% 
of the Yp ro t oso iar, and this number, along with the Y atm0 s 
of Saturn, must be explained in any consistent evolution- 
ary history of the planets. 

Hubbard et al. (1999) showed that the onset of helium 
phase separation in Jupiter leads to model planet ages 
greater than the age of the solar system when Jupiter 
reaches T eff = 124.4 K. However, Hubbard et al. (1999) 
(and Paper I) made the same assumptions outlined in Sec- 
tion 2. They assume all photons that are absorbed are 
done so in adiabatic layers, and that Jupiter's atmospheric 
T — P profile does not deviate from that of an isolated ob- 
ject. From the Galileo entry probe we know that this is 
not the case. In addition, these papers assume Jupiter's 
Bond albedo is a constant 0.343 with time. However, if 
Jupiter's Bond albedo was larger in the past, when H 2 
clouds were higher in the planet's atmosphere, its bond 
albedo could be somewhat larger and less solar radiation 
would be absorbed. This would lead to faster cooling and 
could admit into Jupiter models some phase separation 
before an age of 4.6 Gyr. While evolutionary models con- 
sistently incorporating incident stellar radiation have been 
calculated for HD 209458b (Baraffe et al. 2003; Burrows 
et al. 2003), similar work has not yet been done for Jupiter 
and Saturn. 

4.3. Future Work and Observations 

Since we have yet to work out a consistent evolution- 
ary history for Jupiter and Saturn that explains the at- 
mospheric abundance of helium in both of these planets, 
the possibility that a different element could be undergo- 
ing phase separation must be entertained. This could be 
happening in addition to, or instead of, helium's phase 
separation. Possibilities include neon, which was found in 
only l/10th the solar abundance in Jupiter's atmosphere 
(no value is known for Saturn) and oxygen, since a re- 
cent model of Jupiter's formation (Gautier et al. 2001a, b) 
predicts H 2 may be at least 9.4 times more abundant 



in Jupiter than it is in the Sun. This fairly large abun- 
dance could lead to immiscibility temperatures of several 
thousand K at megabar pressures, but no detailed oxy- 
gen/hydrogen phase diagrams have yet been calculated at 
these pressures. 

If a radial velocity candidate planet is directly imaged in 
the coming years, without an accurate mass determination 
by another method, it will be important to use inhomoge- 
neous evolutionary models that include helium phase sepa- 
ration when trying to deduce the planet's mass from model 
evolution tracks. In addition, if astrometry missions such 
as SIM lead to EGP detections (with mass determinations) 
around stars, it will be worthwhile to include the effects of 
phase separation when trying to understand the sensitivity 
needed to image these planets. A measured luminosity of 
a low-mass EGP with a known mass would be an impor- 
tant observational test of our suggested phase diagram. 
In addition, these larger intrinsic luminosities should be 
taken into account when calculating accurate atmosphere 
models for EGPs, and their corresponding emission and 
reflection spectra. 

While younger planets have the advantage of being sev- 
eral orders of magnitude brighter if their systems are only 
millions of years old, planets found by radial velocity or as- 
trometry perturbations will be attractive targets because 
planetary orbital parameters will already be known. These 
planets that orbit at several AU and beyond around their 
parent stars, if one can use Saturn and Jupiter as a guide, 
may well be ~ twice as luminous as current homogeneous 
models predict. Of the currently available planetary can- 
didates, e Eri b and 55 Cnc d at first glance seem to have 
the greatest likelihood of having helium phase separation 
currently affecting their evolution. However, even though 
e Eri b has a low minimum mass (0.92 Mj (Hatzes et al. 
2000)), its age is likely only 1 Gyr (Drake & Smith 1993). 
55 Cnc b, although it orbits with a semimajor axis of 5.9 
AU, has a minimum mass of 4.05 Mj (Marcy et al. 2002), 
making it too massive to undergo helium phase separation. 

The calculations presented here are an early step in ex- 
ploring phase separation in EGPs. No experimental data 
are yet available for high pressure hydrogen-helium mix- 
tures at planetary temperatures. Advances in understand- 
ing phase separation will occur through high pressure ex- 
periments, theory, further detailed evolutionary models 
of Jupiter and Saturn, and an accurate determination of 
Saturn's atmospheric helium abundance. While we be- 
lieve our calculations for planets less massive than Jupiter 
should be reasonably accurate, since we have calibrated 
the theory to these planets, the extrapolation up to 3.0 
Mj is uncertain, since it is not clear that the immiscibility 
curves we show in Figure 2 maintain this positive slope at 
pressures greater than tens of Mbar. A turnover in the 
immiscibility curves would lead to helium redissolving in 
deeper liquid metallic hydrogen layers, rather than settling 
onto the core. 

5. CONCLUSIONS 

Using the phase diagram described in Paper I, which is 
calibrated to Jupiter and Saturn and allows both planets 
to reach their known ages and effective temperatures, we 
have explored the effects that helium phase separation will 
have on a variety of EGPs. The additional energy liber- 
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ated as helium rains to deeper layers of a planet will signif- 
icantly delay the cooling and contraction of giant planets. 
Once helium phase separation is underway, our inhomoge- 
neous evolutionary models predict luminosities ~ 2 times 
greater than predictions from homogeneous models. This 
will make these giant planets in the 0.15 to 3.0 Mj mass 
range somewhat easier to detect that has been previously 
thought. Improvements in understanding the evolution of 
these EGPs will come through a better understanding of 



our closest giant planets, Jupiter and Saturn. 

Electronic files of the evolutionary models presented 
here can be obtained by contacting JJF. We thank Adam 
Burrows, David Sudarsky, Jonathan Lunine, Ivan Hubeny, 
and Jason Barnes for many interesting conversations while 
this work was underway. JJF is funded by a NASA GSRP 
grant and WBH by NASA PG&G grant NAG5-13775. 
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Fig. 1. — Effective temperature vs. time in Gyr for homogeneous fully adiabatic models of Jupiter and Saturn from Paper I. Jupiter's core 
is 10 Me with a Z eri „ e ; ope =0.059. Saturn's core is 21 Ms with a Z en „ e ; ope =0.030. These values lead to the correct values of each planets' 
radius and moment of inertia at their known effective temperatures. Both planets' known effective temperatures and the age of the solar 
system are shown as dotted lines. Both planets possess Y = 0.27. 
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-1.0 -0.5 0.0 0.5 1.0 1.5 2.0 

log P (Mbar) 

Fig. 2. — Our current understanding of the high pressure phase diagram of hydrogen. Regions of liquid molecular hydrogen (H2) and liquid 
metallic hydrogen (H+) are shown at high temperatures, and their solid counterparts at much lower temperatures. The curve marked PPT? 
is a possible transition from liquid H2 to liquid H + as calculated by Saumon et al. (1995). The solid line marked 50% shows, for the theory of 
Ross (1998), where liquid H2 should be 50% dissociated. Laser shock data points from Collins ct al. (1998) are shown as a dash-dot line. The 
reverberation shock data of Nellis et al. (1999) are shown as pluses (for deuterium) and asterisks (for hydrogen). The large black dot indicates 
the highest pressure that the conductivity of H2 has been measured, which seems to indicate H2 may be 10% dissociated at this point (Nellis 
et al. 1999). The calculated region of helium immiscibility from Hubbard & Dewitt (1985) and Stevenson (1975) is labeled HDW/S. The 
lines marked Y = 0.27 and Y = 0.21 mark the immiscibility boundaries for these two compositions. The parallel lines labeled Pfaff. show 
the helium immiscibility region as calculated by Pfaffcnzeller et al. (1995). Again, the upper boundary is for Y = 0.27 and the lower for 0.21, 
although they are not labeled to avoid clutter. The current internal adiabats of Jupiter and Saturn are shown as heavy lines, while the dashed 
extensions show the pressure range within their cores. Also shown in a heavy line is the adiabat of a hypothetical coreless 0.15 Mj planet 
at an age of 4.5 Gyr. The parallel lines marked PI between the Jupiter and Saturn adiabats is the ad-hoc immiscibility region from Paper 
I. The immiscibility lines are defined to be exactly parallel to the adiabats. The arrow near PI is meant to indicate that the immiscibility 
curves continue at this slope to higher pressures. 
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Fig. 3. — Additional evolutionary models from Paper I including the phase separation of helium from liquid metallic hydrogen. The dotted 
curve for Saturn and the solid curve for Jupiter are the evolutionary tracks shown in Figure 1. The dashed curve includes the phase diagram 
of Hubbard & Dewitt (1985), while the solid curve for Saturn uses the proposed ad-hoc phase diagram from Paper I. The Saturn curves are 
slightly offset at younger ages due to differences in initial core masses and heavy element abundances needed for the each planet to reach 
Saturn's known radius and moment of inertia at 95 K. The Paper I phase diagram allows both Jupiter and Saturn to reach their known ages 
and T e f fS. 
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Fig. 4. — Evolution of the luminosity for planets of mass 0.15 to 3.0 Mj for our standard models with no stellar irradiation and 10 Me 
cores. The dotted lines are models without helium phase separation, while the solid lines include the effects of helium phase separation on 
the planets' cooling. The top curve is for the highest mass planet while the bottom curve is for the lowest mass planet. 
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Fig. 5. — Evolution of T e ff for planets of mass 0.15 to 3.0 Mj for our standard models. The top curve if for the highest mass planet while 
the bottom curve is for the lowest mass planet. 
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Fig. 6. — Evolution of Y a tmos (atmospheric helium mass fraction) for planets of mass 0.15 to 3.0 Mj for our standard models. The lower 
the mass of the planet, the earlier the onset of helium phase separation, and the lower the final Y a tmos at a given age. 
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Fig. 7. — Evolution of the radius of planets of mass 0.15 to 3.0 Mj for our standard models. The top curve is for the highest mass planet 
while the bottom curve is for the lowest mass planet. 
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Fig. 8. — Evolution of the T e ff of irradiated model planets with a constant T e q of 79.2 K. This corresponds to a distance of 10 AU from a 
constant luminosity 1.0 Lq star. Only the lowest mass planets reach their T eq in 10 Gyr. 
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Fig. 9. — The Y a tmos for the model planets at 10 AU. In our theory, when the lowest mass planets reach their T eq , their evolution stalls, 
so no additional helium is lost to deeper layers. The 0.15 and 0.2 Mj planets reach this point. See the text for discussion on this issue. 
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Fig. 10. — Evolution of the T e ff of irradiated model planets with a constant T eq of 112.0 K. This corresponds to a distance of 5 AU from a 
constant luminosity 1.0 Lq star. The evolution of all planets is significantly affected, with the lowest mass planets reaching T eq in less than 
a few Gyr. The 0.15 to 0.5 Mj planets reach their T eq before the onset of helium phase separation. 
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Fig. 11. — The Yatmos for the model planets at 5 AU, up to an age of 20, rather than 10, Gyr. Curve labels are the planets' masses in 
Mj. Due to the stellar heating 0.15 to 0.5 Mj planets reach their T eq before the onset of phase separation, while the 3.0 Mj planet stays 
too warm for 12 Gyr. Once a planet reaches its T eq , its thermal radiation is dominated by thermalizcd stellar photons, rather than intrinsic 
energy. This is most obvious for the 0.7 Mj planet as its Y a tmos nearly reaches an asymptotic value just below 0.24. 
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Fig. 12. — The evolution of isolated planets' radii as a function of time for coreless models and models with heavy element cores of 20 
Mg. Planets are labeled in My. 0.15, 0.2, 0.3, 1.0, and 2.0. For each mass pair, the planet with the larger radius is the coreless model. The 
coreless model planets of mass 0.15 and 0.2 Mj (dashed lines) do not have central pressures high enough for liquid metallic hydrogen to form 
(~ 2 Mbar) at any age, and therefore helium will not become immiscible. 
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Fig. 13. — The evolution of isolated planets' effective temperatures as a function of time for coreless models and models with heavy element 
cores of 20 Mg. Planets are labeled in My. 0.15, 0.2, 0.3, 1.0, and 2.0. For each mass pair, the planet with the lower T e ff is the coreless 
model. The coreless model planets of mass 0.15 and 0.2 Mj (dashed lines) do not have central pressures high enough for liquid metallic 
hydrogen to form (~ 2 Mbar) at any age, and therefore helium will not become immiscible. 
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Fig. 14. — Helium phase separation occurs before or (more likely) during the formation of ammonia clouds, which will delay their settling 
to higher pressures in the planets' atmospheres as they cool. This figure shows, as a function of planet mass and irradiation, how long it 
takes ammonia clouds to reach 1 bar pressure in the planets' atmospheres. If the effects of phase separation are included, planets at 5 AU 
will never form ammonia clouds that reach 1 bar. 



